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Device independent dimension witnesses provide a lower bound on the dimensionality of classical 
and quantum systems in a "black box" scenario where only correlations between preparations, 
measurements and outcomes are considered. We address the problem of the robustness of dimension 
witnesses, namely that to witness the dimension of a system or to discriminate between its quantum 
or classical nature, even in the presence of loss. We consider the case when shared randomness 
is allowed between preparations and measurements and we provide a threshold in the detection 
efficiency such that dimension witnessing can still be performed. 



I. INTRODUCTION 

For several experimental setups, a description that 
completely specifies the nature of each device is unsatis- 
factory. For example, in a realistic scenario the assump- 
tion that the provider of the devices is fully reliable is 
often overoptimistic: losses unavoidably affect the imple- 
mentation, thus turning it away from its ideal descrip- 
tion. A device independent description of an experimen- 
tal setup does not make any assumption on the involved 
devices, which are regarded as "black boxes" , while only 
the knowledge of the correlations between preparations, 
measurements and outcomes is considered. In this sce- 
nario, a natural question is whether it is possible to derive 
some properties of the non-characterized devices instead 
of assuming them, building only upon the knowledge of 
these correlations. In general one could be interested in 
bounding the dimension of the systems prepared by a 
non-characterized device; one could also ask whether a 
source is intrinsically quantum or can be described clas- 
sically. The framework of device independent dimension 
witnesses (DIDWs) provides an effective answer to these 
questions, suitable for experimental implementation and 
for application in different contexts, such as quantum 
key distribution (QKD) or quantum random access codes 
(QRACs). 

DIDWs were first introduced in [l[ in the context of 
non-local correlations for multi-partite systems. Subse- 
quently, the problem of DIDWs was related to that of 
QRACs in [2[, and in [3[ it was reformulated from a dy- 
namical viewpoint allowing one to obtain lower bounds 
on the dimensionality of the system from the evolution 
of expectation values. A general formalism for tackling 
the problem of DIDWs in a prepare and measure scenario 
was recently developed in [J]. The derived formalism al- 
lows one to establish lower bounds on the classical and 
quantum dimension necessary to reproduce the observed 
correlations. Shortly after, the photon experimental im- 
plementations followed, making use of polarization and 
orbital angular momentum degrees of freedom [5] or po- 
larization and spatial modes |6| to generate ensembles of 
classical and quantum states, and certifying their dimen- 
sionality as well as their quantum nature. 



DIDWs also allow reformulating several applications in 
a device- independent framework. For example, dimen- 
sion witnesses can be used to share a secret key between 
two honest parties. In Q, the authors present a QKD 
protocol whose security against individual attacks in a 
semi-device independent scenario is based on DIDWs. 
The scenario is called semi-device independent because 
no assumption is made on the devices used by the honest 
parties, except that they prepare and measure systems 
of a given dimension. Another application is given by 
QRACs, that make it possible to encode a sequence of 
qubits in a shorter one in such a way that the receiver 
of the message can guess any of the original qubits with 
maximum probability of success. In [1, @ QRACs were 
considered in the semi-device independent scenario, with 
a view to their application in randomness expansion pro- 
tocols. 

Clearly any experimental implementation of DIDWs 
is unavoidably affected by losses - that can be modeled 
as a constraint on the measurements - and can reduce 
the value of the dimension witness, thus making it im- 
possible to witness the dimension of a system. These 
considerations make of the utmost importance the ques- 
tion whether it is possible to perform reliable dimension 
witnessing in realistic scenarios and, in particular, with 
non-optimal detection efficiency. We refer to this prob- 
lem as the robustness of device independent dimension 
witnesses. Despite its relevance for experimental imple- 
mentations and practical applications, this problem has 
not been addressed in previous literature. The aim of this 
work is to fill this gap. We consider the case where shared 
randomness between preparations and measurements is 
allowed. Our main result is to provide the threshold in 
the detection efficiency that can be tolerated in dimen- 
sion witnessing, in the case where one is interested in 
the dimension of the system as well as in the case where 
one's concern is to discriminate between its quantum or 
classical nature. 

The paper is structured as follows. In Section |TT] we 
discuss some relevant properties of the sets of quantum 
and classical correlations and introduce the concept of 
dimension witness as a tool to discriminate whether a 
given correlation matrix belongs to these sets. In Section 
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IIIII we provide a threshold in the detection efficiency that 
is allowed in witnessing the dimensionality of a system 
or in discriminating between its classical or quantum na- 
ture, as a function of the dimension of the system. We 
summarize our results and discuss some further develop- 
ments - such as dimension witnessing in the absence of 
correlations between preparations and measurements or 
entangled assisted dimension witnessing - in Sect. IIVI 



II. 



DEVICE INDEPENDENT DIMENSION 
WITNESSES 



Let us first fix the notation [10j . Given an Hilbert 
space %, we denote with LinH the space of linear oper- 
ators X : % — > %. A quantum state in % is represented 
by a density matrix, namely a positive semi-definite ma- 
trix p £ lAnH such that Tr[/o] = 1. Given a pure state 
\ip) £ H, we denote with ip := \ip)(ip\ the corresponding 
projector. A set R = {p^} of states is said to be classical 
when the states commute pairwise, [pi,Pk] = for any 
i, k. A general quantum measurement is represented by 
a Positive-Operator Valued Measure (POVM), namely a 
set of positive semi-definite Hermitian matrices IP such 
that £\ IP = I. A POVM n = {IP} is said to be clas- 
sical when [IP, IF] = for any j, I. The joint probability 
of outcome j given input state pi is given by the Born 
rule, namely Pj\i = Tr[piU j ]. 

The general setup introduced in Q for performing 
device independent dimension witnessing is given by a 
preparing device (let us say on Alice's side) and a mea- 
suring device (on Bob's side) as in Fig. [T] In the most 
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FIG. 1. (Color online) Setup for witnessing the dimension of 
a quantum or classical system. In the most general scenario 
considered here, Alice and Bob share a hidden random vari- 
able A. Alice (on the left hand side) owns a preparing device 
which sends the state p it \ to Bob whenever Alice presses but- 
ton i £ [1, M]. Bob owns a measuring device that performs 
measurement II^a on the received state whenever Bob presses 
button k £ [1, K], giving the outcome j £ [1, AT]. 

general scenario, the devices may share a priori correlated 
information, classical and quantum. However, in many 
realistic situations, one can assume that the preparing 
and measuring devices are uncorrelated and that all the 
correlations observed between the preparation and the 
measurement are due to the mediating particle connect- 
ing the two devices. An intermediate and also valid pos- 
sibility is to assume that the devices only share classical 



correlations. In this case, the value of a random variable 
A distributed according to q\ is accessible to preparing 
and measuring devices. In this work we focus on this last 
possibility. Alice chooses the value of index i € [1, M] and 
sends a fixed state pi t \ £ LinH to Bob. Bob chooses the 
value of index k £ [l,K] and performs a fixed POVM 
Hk,\ on the received state, obtaining outcome j £ [1, AT]. 
After repeating the experiment several times (we con- 
sider here the asymptotic case) , they collect the statistics 
about indexes i,j,k obtaining the conditional probabili- 
ties PjUfe. Note that we also implicitly assume that we 
are dealing with independent and identically distributed 
events. 

We now introduce the set Q (the set C) of correlations 
achievable with quantum (classical) preparations. 

Definition 1 (Set of quantum correlations). For any 
M, K, V, d £ N we define the set of quantum correlations 
Q(M,K,N,d) as the set of correlations Vj\i,k with i £ 
[1,-M], k £ and j £ [1,AT] such that there exist 

a Hilbert space % with dim"H = d, a quantum set R = 
{pi £ LinH}f of states and a set P = {ILjf of POVMs 
Hk = {IIj. £ L'mH}^ for which PjU t k = Tr^IL^,], namely 

Q := {p | 3 d- dimensional Hilbert space %, 

3 quantum set {pi £ LinH}^ of states, 
3 set {n fe }f of POVMs II fe = {n{ £ LinHl^ 

Pj\i,k 



N 



such that h — Tr[piLT|,]}. 



Definition 2 (Set of classical correlations). For any 
M, K,N,d £ N we define the set of classical correla- 
tions C(M, K, N,d) as the set of correlations Pj\i_k with 
i £ [1, M], k £ [lj -K] and j £ [1, N] such that there exist 
a Hilbert space % with dim"H = d, a classical set R = 
{pi £ UxiU}¥ of states and a set P = {ILjf of POVMs 
lift = {LT^, £ LinH}^ for which pj\ i k = Tr^LI^], namely 

C := {p \ 3 d-dimensional Hilbert space H, 

3 classical set {pi £ \AnT-C\i of states, 

3 set {n fc }f of POVMs II fe = {n J fe £ LinH}f 

such that Pj\i,k — Tr[pill^]}. 

We write Q and C omitting the parameters M, K, N, d 
whenever they are clear from the context. 

Remark 1. We notice that, when shared randomness 
is allowed between quantum (classical) preparations and 
measurements, the set of achievable correlations is qiven 
by Conv Q (ConvC), where for any set X we denote with 
Coxrv X the convex hull of X . 

The following Lemma shows that it is not restrictive 
to consider only classical POVMs in the definitions of 
classical correlations. 

Lemma 1. For any correlation p = {pj\i.k} G C there 
exist a classical set R = {pi} of states and a set Q = 
{Afe} of classical POVMs Ak = {Aj.} such that PjU t k = 
Tr[pj Aj <k ] . 
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Proof. By hypothesis there exist a classical set R = {pi} 
of states and a set P = {II fe } of POVMs U k = {u{} 
such that PjUk — TrL^n^] for any i,j,k. Take Aj, = 
where {\i}} is an orthonormal basis with 
respect to which the pi's are diagonal (it is straightfor- 
ward to verify that A k > for any k, j and JV Aj, = 7 for 

any fc). We have Vj\i,k = Tr[piA^] for any i,j,k, which 
proves the statement. □ 

Since classical correlations can always be reproduced 
by quantum ones, we immediately have C C Q and 
ConvC C Conv Q. Moreover, by definition we have 
C C ConvC and <2 C Conv Q. Here we show an example 
where C is non-convex (namely C C ConvC) and C C Q. 
Take M = 3, K = 2, AT = 2, and d = 2 and consider the 
following conditional probability distribution 



pm 




Pi|i,2 




(1) 



First we show that j? € ConvC. Indeed p can be 
obtained when Alice and Bob share classical correla- 
tions represented by a uniformly distributed random vari- 
able A taking values 1 , 2 making use of the classical set 
R = {Pi,x} of states and of the set P = {Hk, x} of classical 
POVMs n fciA = {U{ J, with 



Pi, 2 



P2,l 
P2,2 



|1>(1| 



P3,l 
P3,2 



and 



ni,i = 
nf a = 



n^ = |o)<o|, 
nL = |i)<i|, 



which proves that p e ConvC. 

Now we show that p£ Q. Indeed p can be obtained by 
Alice and Bob making use of the quantum set R = {/?$} 
of states and of the set P = {Hk} of quantum POVMs 
U k = {IPJ, with 



pi = |o)(o|, p 2 = |+)(+|, P3 = |i)(i|, 



and 



ni = |o)<o|, nl = \+)(+\, 

which proves that p G Q. 

Finally, we verify that if Alice and Bob make use of 
classical sets of states and POVMS and do not have ac- 
cess to shared randomness there is no way to achieve 
the probability distribution p given by Eq. ([1]). Indeed, 
to have perfect discrimination between p\ and ps with 
POVM III [see Eq. JIJ], one must take p\ and p$ orthog- 
onal - let us say without loss of generality p\ = |0)(0| 
and p 3 = |1)(1|, and n} = |0)(0| and II? = |1)(1|. Due 
to the hypothesis of classicality of the sets of states, p2 



must be a convex combination of p\ and p 3 . Then, in 
order to have Pj\2,i as in Eq. ([1]), one has to choose 
92 = (pi + P3)/2 = 1/2. Finally, the only possible choice 
for H2 is IIj = I and = 0, which is incompatible with 
the remaining entries of Pju% in Eq. ([1]). This proves 
that p g C. 

The relations between the sets of quantum and classical 
correlations are schematically depicted in Fig. [2j 




FIG. 2. (Color online) Schematic representation of the sets 
of classical and quantum correlations between preparations, 
measurements and outcomes. Dashed line represents the 
(non-convex) set C of classical correlations without shared 
randomness; the rectangle represents the set ConvC of clas- 
sical correlations with shared randomness; the ellipsoid rep- 
resents the set Conv Q of quantum correlations with shared 
randomness. 

We can now introduce DIDWs. Building only on the 
knowledge of Pj\i t k, our task is to provide a lower bound 
on the dimension d of T~L or to certify that the states pi 
have to be quantum if their dimension is assumed to be 
smaller than a given value. 

Definition 3. For any set of correlations X between M 
preparations and K measurements with N outcomes, a 
device independent dimension witness Wx{p) is a func- 
tion of the conditional probability distribution p — {pj\i t k} 
with i S [1, M], fee [1 , iV] 7 and j € [1, N] such that 



W x (p)> L^p^X, 
for some L which depends on Wx- 



(2) 



Motivated by Remark [Q for any M,N,K,d £ 
N when X = Conv C(M, N, K, d) [when X = 
Conv Q(M, N, K,d)} we say that Wx(p) is a classical 
(quantum) dimension witness for dimension d in the pres- 
ence of shared randomness. Given a set R = {pi,x} 
of states and a set P = {H- k ,\\ of POVMs Il k ^ = 
{n{ <x }, we define Wc nvc(R, P) ■= Wconv c(p) with 
P = {pj\i,k} and pj\ i>k = J^x 9a Traill], and analogously 

for WconvQ- 

In this work we will consider only linear DIDWs, 
namely inequalities of the form of Eq. ([2J such that 



W(p) := c-p= ^2ci^kPj\i,k> 



(3) 



where c is a constant vector. 
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Notice that for any function W(p) and constant L, the 
witness W(p) > L is only a representative of a class of 
equivalent witnesses such that if W'(p) > V is a member 
of the class, then W(p) > L if and only if W'(p) > V 
for any conditional distribution p. The following Lemma 
provides a transformation that preserves this equivalence. 

Lemma 2. Given a function W{p) = ■ k Ci,j,kPj\i,k 
and a constant L, take W'(p) = j k c i j kPj\i,k with 
= c ijjtk + a it k and L' = L + £\ fe a i>k for any 
cti^k that does not depend on outcome j. Then one has 
W(p) > L if and only if W'(p) > L' for any p. 

Proof. It follows immediately by direct computation. □ 

In the following our task will be to find a set R of 
quantum states and a set P of POVMs such that a linear 
witness W(R,P) maximally violates inequality ([2]). The 
following Lemma allows us to simplify the optimization 
problem. 

Lemma 3. The maximum of any linear dimension wit- 
ness W(R, P) is achieved by an ensemble R of pure states 
and without shared randomness. 

Proof. The thesis follows immediately from linearity. □ 

Due to Lemma[3jthe maximization of Eq. ([3]) is equiv- 
alent to the maximization of 

w{r,p) = Y;^ 3 MK\^)i 

i,j,k 

over the sets R = {tpi} of pure states and the sets P — 
{U k } of POVMs n fc = {!!{}. 



III. ROBUSTNESS OF DIMENSION 
WITNESSES 

In practical applications, losses (due to imperfections 
in the experimental implementations or artificially intro- 
duced by a malicious provider) can noticeably affect the 
effectiveness of dimension witnessing. The main result of 
this Section is to provide a threshold value for the de- 
tection efficiency which allows to witness the dimension 
of the systems prepared by a source or to discriminate 
between its quantum or classical nature, when shared 
randomness between preparing and measuring devices is 
allowed. 

The task is to determine whether a given conditional 
probability distribution belongs to a particular convex 
set, namely ConvC or Conv Q (see Remark Q]). The sit- 
uation is illustrated in Figure [3] The experimental im- 
plementation is constrained to be lossy, namely it can be 
modeled considering an ideal preparing device followed 
by a measurement device with non-ideal detection effi- 
ciency. This means that any POVM Hk,\ on Bob's side 




FIG. 3. (Color online) The Figure illustrates the problem of 
the robustness of device independent dimension witness. The 
convex hulls Conv Q and Conv C of the sets of quantum and 
classical correlations are represented as in Figure [2] In the 
presence of loss, only a subset of the possible correlations is 
attainable. The subset, surrounded by bold line in the figure, 
is parametrized by detection efficiency n. The task is to find 
the threshold value in n such that dimension witnessing is 
still possible. For example, when the task is to discriminate 
between the quantum or classical nature of a source, one is 
interested in achieving correlations in the dashed region, and 
our goal is to determine the values of n such that this area is 
not null. 

is replaced by a POVM Tii \ with detection efficiency rj, 
namely 

n^:={„n M ,(i-^)j}. (4) 

We notice that each lossy POVM has one outcome more 
than the ideal one, corresponding to the no-click event. 
In a general model, the detection efficiency t] may be 
different for any POVM Hk.\- Nevertheless, in the fol- 
lowing we assume that they have the same detection ef- 
ficiency, which is a reasonable assumption if the detec- 
tors have the same physical implementation [111 ]. Anal- 
ogously given a set P = {Ilk, a} of POVMs we will de- 
note with = {ni jl} the corresponding set of lossy 
POVMs. Upon defining pW := {p^ J with p ( ^ i k = 
J2\ q\ ^\Pi,)J^k,\\ one clearly has 

pW=r,pW+(l-)))p to l (5) 

To attain our task we maximize a given dimension wit- 
ness over the set of lossy POVMs as given by Eq. Q. 
Due to the model of loss introduced in Eq. (H|) and to 
the freedom in the normalization of dimension witnesses 
given by Lemma [21 in the following without loss of gen- 
erality for any dimension witness W as given in Eq. ([3]) 
it is convenient to take 

Ci,N,k = 0, Vi,fc. (6) 

Then we have the following Lemma. 

Lemma 4. Given a set R = {pi,\]iL\ of states and a set 
P = {n fe , A }f =1 of POVMs n k ,x = {IL 7 Jf^ 1 , for any 
linear dimension witness W(p) — ^ k Ci,j,kPj\i,k with 
i G [1, Af], j G [X, JV], and k G normalized as in 

Eq. §6§ one has 

W(R,P M ) =r 1 W(R,P ( - 1) ). 
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Proof. One has 

W(R, PW) = J2 Mm + (1 - V)pf\l k 
= nW(R,P {1) ), 



where the first equality follows from Eq. §5§ and the 
second from the fact that W(p^) = due to the nor- 
malization given in Eq. ([6]). □ 

In particular from Lemma 0] it follows that for any 
linear dimension witness W one has 

max P^) = rimaxWiKP^), 

R,P B„P 

arg max W(R, P^) 

R 7 P R.P 

Due to Lemma HI it is possible to recast the optimiza- 
tion of dimension witnesses in the presence of loss to the 
optimization in the ideal case. Then due to Lemma [3] 
it is not restrictive to carry out the optimization with 
pure states and no shared randomness. Consider the case 
where M — d + 1, K = d, and N = 3. Using the tech- 
nique discussed in Appendix [X] one can verify that the 
witness given by Eq. Q with the following coefficients 



-1 iii + k<M,j = l 
+1 if i + k = M + l, j 
otherwise 



(7) 



is the most robust to non-ideal detection efficiency. This 
fact should not be surprising, as we notice that this wit- 
ness relies on only 2 out of 3 outcomes. According to 
we denote it Id+i- In [U (see also [l3j]) it was conjectured 
that for any dimension d the dimension witness I c i+i is 
tight in the absence of loss. 

Now we provide upper and lower bounds for the max- 
imal value I^+i '■= max fl,p Id+i where the maximiza- 
tion is over any set R — {pi G Lin"H} of states and 
any set P = {n fc } of POVMs II fe = {H j k G Lin"H} with 
dim'H = d. 

Lemma 5. For any dimension d we have > + 

Proof. The statement follows from the recursive expres- 
sion Id+i = I<i + C, where 



c -.= -J^iiniivo + (^+i|n^+i) 



and noticing that Id and C can be optimized indepen- 
dently. □ 

A tight upper bound for I 3 was provided in Q . In the 
following Lemma we provide a constructive proof suitable 
for generalization to higher dimensions. 



Proof. The statement follows from standard optimiza- 
tion with Lagrange multipliers method and from the 
straightforward observation that given two normalized 
pure states \vq) and \vi), if a pure state \u) can be de- 
composed as follows 



\u) = (vq\u)\vq) + (vi\u)\Vi), 

then \{vq\u)\ = \(vi\u)\. 



□ 



Making use of Lemmas [5] and El we provide upper and 
i^.j as follows 

d - 2 + V2 < I* d+1 < d, (8) 



lower bounds on as follows 



where the second inequality follows from the non discrim- 
inability of d + 1 states in dimension d (see [J]). 

We now make use of these facts to provide our main re- 
sult, namely a lower threshold for the detection efficiency 
required to reliably dimension witnessing. We consider 
the problem of lower bounding the dimension of a system 
prepared by a non-characterized source in Proposition 
[TJ as well as the problem of discriminating between the 
quantum or classical nature of a source in Proposition [21 

Proposition 1. For any d there exists a dimension 
witnessing setup such that it is possible to discrimi- 
nate between the quantum and classical nature of a d- 
dimensional system using POVMs with detection effi- 
ciency rj whenever 



V > Vgc ■= (d - 1)/Jd+1 



Furthermore one has 
d-1 



< 



Vqc 



< 



d-1 



d-2 + V2' 



(9) 



(10) 



Proof. We provide a constructive proof of the statement. 
Take M = d + 1, K = d, and N = 3, and we show that 
Id+i satisfies the thesis. 

We notice that the maximum value of Id+i attainable 
with classical states is given by d — 1 [|J . Then ?7 qc is the 
minimum value of the detection efficiency such that Id+i 
can discriminate a quantum system from a classical one. 

Due to Lemma [4] we have Eq. ([9]) . From Eq. ([§]) 
the lower and upper bounds for n qc given in Eq. (|10[) 
straightforwardly follow. □ 

Notice that Id+i in Eq. © can be numerically evalu- 
ated with the techniques discussed in Appendix [Bj Fig- 
ure |H plots the value of ry qc for different values of the 
dimension d of the Hilbert space T-L. The threshold in 
the detection efficiency when d — 2 is n qc = l/v2 ; going 
asymptotically to 1 with d as ~ 1 + 1 / d. 

Proposition 2. For any d there exists a dimension wit- 
nessing setup such that it is possible to lower bound the 
dimension of a d + 1-dimensional system using POVMs 
with detection efficiency rj whenever 



Lemma 6. For dimension d = 2 we have /■? = v2- 



h+i/d. 



(11) 
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FIG. 4. (Color on line) Threshold value (middle line) of the 
detection efficiency r] qc as in Eq. D as a function of the 
dimension d, obtained through numerical optimization of Id+i 
with Algorithm [2] The lower bound (lower line) and upper 
bound (upper line) given by Eq. (|l(jp are also plotted. As 
expected, the upper bound is tight for d = 2. The detection 
efficiency -q qc asymptotically goes to 1 as d —¥ oo since its 
upper and lower bound do the same. 



FIG. 5. (Color on line) Threshold value (upper line) of the 
detection efficiency jjdim as in Eq. (j 1 1 p as a function of the 
dimension d, obtained through numerical optimization of Id+i 
with Algorithm [2] The lower bound (lower line) given by Eq. 
(|12p is also plotted. As expected, the lower bound is tight for 
d = 2. The detection efficiency ?7dim asymptotically goes to 1 
as d — s> oo since its lower bound does the same (and Tjdim < 1 
is a trivial upper bound). 



Furthermore one has 

2 \/2 

ridim > 1 - , ■ (12) 
a 

Proof. We provide a constructive proof of the statement. 
Take M = d + 1, K = d, and N = 3, and we show that 
Id+i satisfies the thesis. 

We notice that the maximum value of Id+i attainable 
in any dimension > d is given by d Then rydim is the 
minimum value of the detection efficiency such that Id+i 
can lower bound the dimension of a d + 1 dimensional 
system. 

Due to Lemma 0] we have Eq. (fPTj) . From Eq. © the 
lower bound to rydim given by Eq. (fT2"j) straightforwardly 
follows. □ 

Notice that Id+i in Eq. (jlll) can be numerically eval- 
uated with the techniques discussed in Appendix [SJ Fig- 
ure O plots the value of ?7dim for different values of the 
dimension d of the Hilbert space W. The threshold in 
the detection efficiency when d = 2 is ?; qc = 1/V2, going 
asymptotically to 1 with d as ~ 1 + 1/d. We notice that 
%im grows faster than r) qc , thus showing that for fixed 
dimension, the discrimination between the quantum or 
classical nature of the source is more robust to loss than 
lower bounding the dimension of the prepared states. 

IV. CONCLUSION 

In this work we addressed the problem whether a lossy 
setup can provide a reliable lower bound on the dimen- 



sion of a classical or quantum system. First we pro- 
vided some relevant properties of the sets of classical and 
quantum correlations attainable in a dimension witness- 
ing setup. Then we introduced analytical and numerical 
tools to address the problem of the robustness of DIDWs, 
and we provided the amount of loss that can be tolerated 
in dimension witnessing. The presented results are of the 
utmost importance for experimental implementations of 
DIDWs, and can be naturally applied to semi-device in- 
dependent QKD and QRACs. 

We notice that, while we provided analytical proofs of 
our main results, i.e. Propositions Q] and [21 their optimal- 
ity as a bound relies on numerical evidences. In particu- 
lar, they are optimal if the dimension witness Id+i is in- 
deed the most robust to loss for any d, which is suggested 
by numerical evidence obtained with the techniques of 
Appendix [S] and Appendix [Bj Thus, a legitimate ques- 
tion is whether the bounds provided in Propositions Q] 
and[5]are indeed optimal. Moreover, it is possible to con- 
sider models of loss more general than the one considered 
here, e.g. one in which a different detection efficiency is 
associated to any POVM. 

A natural generalization of the problem of DIDWs, in 
the ideal as well as in the lossy scenario, is that in the 
absence of correlations between the preparations and the 
measurements. In this discussed in this work, 

the relevant sets of correlations are Q and C, which are 
non-convex as shown in Section [TTJ The non convexity 
of the relevant sets allows the exploitation of non-linear 
witnesses - as opposed to what we did in the present 
work. An intriguing but still open question is whether 
there are situations in which this exploitation allows to 
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dimension witness for any non-null value of the detection 
efficiency. 

Another natural generalization of the problem of 
DIDWs is that of entangled assisted DIDWs, namely 
when entanglement is allowed to be shared between the 
preparing device on Alice's side and the measuring device 
on Bob's side. This problem is similar to that of super- 
dense coding [I4| . Consider again Fig. [T] In the simplest 
super-dense coding scenario, Alice presses one button out 
of M = 4, while Bob always performs the same POVM 
[K = 1) obtaining one out of N = 4 outcomes. The 
dimension of the Hilbert space W is dim('H) = 2, but 
a pair of maximally entangled qubits is shared between 
the parties. In this case, the results of [l4j imply that a 
classical system of dimension 4 (quart) can be sent from 
Alice to Bob by sending a qubit (corresponding to half 
of the entangled pair). 

Consider the general scenario where now the two 
parties are allowed to share entangled particles. The 
super-dense coding protocol automatically ensures that 
by sending a qubit Alice and Bob can always achieve 
the same value of any DIDW as attained by a classi- 
cal quart. Remarkably, the super-dense coding proto- 
col turns out not to be optimal, as we identified more 
complex protocols beating it. In particular, we found 
a (M = 4, K — 2, N — 4) situation for which, upon 
performing unitary operations on her part of the entan- 
gled pair and subsequently sending it to Bob, Alice can 
achieve correlations that can not be reproduced upon 
sending a quart. This thus proves the existence of com- 
munication contexts in which sending half of a maximally 
entangled pair is a more powerful resource than a classi- 
cal quart. This observation is analogous to that done in 
@, where it was shown that entangled assisted QRACs 
(where an entangled pair of qubits is shared between the 
parties) outperform the best of known QRACs. For these 
reasons we believe that the problem of entangled assisted 
DIDWs deserves further investigation. 



Algorithm 1. For any set R^> = {ip\ } of pure states 
and any set P<°) = {flf'} of POVMs U k 0) = {n£ (0) }, 



1. let 



(l-e)/ + e E Jlfc cwX ,(B) ' W), 
(l-e)/ + eEiCij lfc ^ (n) 



2. let n 



3,(n+l) 



n 



j>(«) 



3. normalize \ip. 



(n+1) \ = ll,7, (n+1) ll- 1 /2|^(™+ 1 ) 



normalize H 



E 



j,(n+l) 



with S k 



As for all steepest-ascent algorithm, there is no protec- 
tion against the possibility of convergence toward a local, 
rather than a global, maximum. Hence one should run 
the algorithm for different initial ensembles in order to 
get some confidence that the observed maximum is the 
global maximum (although this can never be guaranteed 
with certainty). Any initial set of states and any initial 
set of POVMs can be used as a starting point, except 
for a subset corresponding to minima of W(R, P). These 
minima are unstable fix-points of the iteration, so even 
small perturbations let the iteration converge to some 
maxima. The parameter e controls the length of each it- 
erative step, so for e too large, an overshooting can occur. 
This can be kept under control by evaluating W(R, P) at 
the end of each step: if it decreases instead of increasing, 
we are warned that we have taken e too large. 

Referring to Fig. [TJ the simplest non-trivial scenario 
one can consider is the one with M = 3 preparations 
and K = 2 POVMs each with N = 3 outcomes, one of 
which corresponding to no-click event. In this case one 
has several tight classical DIDWs. Applying Algorithm 
[T]we verified that among them the most robust to loss is 
given by Eq. ^ with coefficients given by Eq. ([7}. 
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Appendix A: Numerical optimization of dimension 
witnesses 

Given a linear dimension witness W the following al- 
gorithm converges to a local maximum of W(R, P). 



Appendix B: Numerical optimization of Id+i 

The following Lemma proves that the POVMs maxi- 
mizing Id+i for any dimension d are such that one of their 
elements is a proj ector on a pure state, thus generalizing 
a result from (l5j . 

Lemma 7. For any dimension d, the maximum of Id+i 
is achieved by a set P = {U k } of POVMs U k = {W k } 
with a projector with ranknj. = 1 for any k. 

Proof. For any fixed set R = {ipi} of pure states de- 
fine A k := -J^i^k^' B : = i>k, and x k ■= M + B k . 
Then clearly A k < 0, B k > and rankS^ = 1 for any 
k. From Eq. ([3]) it follows immediately that the op- 
timal set P* = {n*} of POVMs n* = {U* k j } is such 
that n^ 1 = argminjji Tr[XII^]. The optimum of Ij+i is 
achieved when is the sum of the eigenvectors of X k 
corresponding to positive eigenvalues. 
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Upon denoting with Xi(Ak) > ■■■ > \ n {Ak) the eigen- 
values of Ak, the Weyl inequality (see for example jig ]) 
Ai(Xfc) < \i(Ak) + \ n (Bk) holds for any n. Since 
^i(Ak) < and \ n {B}f) — for any k and for any n^O, 
the thesis follows immediately. □ 

Algorithm[T]can be simplified using Lemma[7l The fol- 
lowing algorithm converges to a local maximum of Id+i- 

Algorithm 2. For any set R^> = {ip^} of pure states 
and any set p(°) = {nj; 0) } of POVMs = {Il£ (0) }, 

1. let\^ +1) ) = l^+eEi.fcCy.fc^l^)^), 

2. hi^^i^^ + .e^c^^V^)!^), 



3. normalize |V|" +1) ) = ||# +1) ||" 1/2 |# +1) ), 
I normalize \n { k n+1) ) = | !vr^ +1) 1 1 ~ X/2 |7f£" +1) > . 

The same remarks made about Algorithm [T] hold true 
for Algorithm [2] Nevertheless, we verified that in prac- 
tical applications Algorithm [2] always seems to converge 
to a global, not a local maximum. This can be explained 
considering that without loss of generality it optimizes 
over a smaller set of POVMs when compared to Algo- 
rithm [U Moreover, we noticed that the optimal sets of 
states and POVMs are real, namely there exists a basis 
with respect to which states and POVM elements have 
all real matrix entries. A similar observation was done in 
[jj} in the context of Bell's inequalities. 
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